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ABSTRACT 

The  occurrence  of  open-ocean  deep  convection  require^  a  background  cyclonic  circu¬ 
lation  and  a  preconditioning.  Both  Ci»nditions  reduce  the  siiatincalion  t>r  the  water  col¬ 
umn  within  the  cyclonic  gyre,  which  will  then  become  eligible  for  convection  if  the  .surface 
forcing  is  sulTicicntly  intense.  Therefore,  for  open-tKcan  <lecp  cttnvcciion  the  gcncr.ition 
of  vertical  cells  inside  the  vortex  b\'  any  mechanism  (not  necessary  bv  pure 
ihcrmotlynamical  prttccsscs')  is  ultimately  important.  There  .arc  tyvo  dtnamicai  mech¬ 
anisms  for  inducing  vcriic.al  cells  inside  a  stably  siraiifietl  vonc.x:  baroclinic  and 
centrifugal  instabilities.  The  combination  of  the  two  is  called  symmetric  instability.  In 
order  to  verify  the  imprvrtancc  of  symmetric  in.stabiliiy  on  the  generation  of  vertical  cells 
inside  the  vortex,  a  tangential  velocity  Held  with  Gaussian  slisiribution  in  both  radial  and 
vertical  under  stable  stratification  is  taken  as  a  mean  flow  field.  The  disturbances 
produced  from  this  mean  Row  are  assumed  to  be  two  dimensional  (no  azimuthal  de¬ 
pendency)  and  described  by  a  strcamfunciion  in  the  radial-vcmcal  sections.  This 
strcamfunciion  satisfies  a  second-order  partial  differential  equation.  The  numerical  sol¬ 
utions  show  the  generation  of  vertical  cells  inside  the  vortex.  The  strength  and  structure 
of  these  cells  largely  depend  on  the  four  parameters:  Burger  number  Br/  = 

Rossby  number  /?«  =  baroir<>pic  index  m,.  and  baroclinic  index  rn..  The  larger  the 
values  of  Bo.  m„  and  ,  or  the  smaller  the  value  of  Bu  (weaker  stratification  for  a  given 
size  of  vortex),  the  stronger  the  vertical  circulation.  The  time  r.atc  change  of  density 
(density  redistribution)  gcncratcti  by  a  vortex  and  horizontally  .ivcragctl  inside  the  vortex, 
indicates  the  decrease  of  density  in  the  lower  part  of  the  vortex,  and  the  increase  of  den¬ 
sity  in  the  upper  part  of  the  vortex.  This  prfrccss.  on  the  time  scale  of  weeks,  decreases 
the  static  st.Tbility  and  serves  ,qs  one  contributor  to  prcconditiiming  for  the  open  ocean 
deep  convection. 


INTRODUCTION 

The  association  of  deep  convcctiim  with  strong  vortices  has  long  been  observed. 
Gordon  (I97S)  reported  the  occurrence  rT  deep  convection  in  February  1977  within  the 
central  region  of  the  Wcriricll  Gyre  west  of  Maud  Rise,  A  column  of  w-atcr  (14  km  radius) 
was  ob.scrvcd  in  which  the  normal  Antarctic  .stratification  of  temperature  was  absent.  The 
column  appeared  a.s  a  warm,  low-salinity  surface  layer  occupying  approximately  the  up¬ 
per  190  m  above  a  cold  nearly  homogeneous  w'aicr  column,  with  a  cyclonic  eddy  (surface 
velocity  greater  than  .'iO  cm.s)  extending  to  at  least  4000  m. 
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One  may  ask  such  questions:  is  there  any  relation  bcueccn  the  vortex  and  the  deep 
convection?  and  what  is  the  relation?  The  main  task  of  this  paper  is  to  answer  these 
questions  and  to  present  a  possible  mechanism  for  generating  vertical  cells  by  a  vortex. 
In  a  gcncr.ll  sense,  convection  can  be  defined  as  vertical  overturning  (Chu.  IW|),  no 
matter  whether  the  process  that  generated  it  is  thermodynamical  or  dynamical. 

Vortices  represent  the  most  energetic  components  of  the  mcso-scnic  eddy  field  in  the 
world  oceans.  Their  large  azimuthal  swirl  speeds,  compared  with  the  rate  at  which  they 
translate,  allows  them  to  carry  substantial  volumes  of  fluid  with  them  (ivcr  long  periods 
of  time.  Rayleigh  (ISXO,  |9l6l  first  investigated  the  instability  of  rotating  structures  in 
non*rotating  fluids.  He  conskicrcti  a  basic  swirling  flow  of  an  inviscUl  fluid  which  moves 
with  angular  vcIcKity  an  arbitrary  function  of  the  distance  r  from  the  axis  of  rota¬ 
tion.  By  a  simple  physical  argument  Rayleigh  then  derived  his  celebrated  criterion  for 
stability.  Rayleigh  s  circulation  criterion  states  that  a  necessary  and  sufficient  condition 
for  stability  to  axisymmctric  slisturbances  is  that  the  square  of  the  angular  momentum 
docs  not  decrease  radially  anywhere. 


1  d 


(I) 


where  the  function,  C(r),  Is  called  the  Rayleigh  Discriminant.  The  vortex  can  be  either 
ccnirifugally  stable  or  unstable,  ilcpending  upon  the  sign  of  the  Rayleigh  Discriminant. 
C.  Rayleigh  further  observed  that  there  Is  an  analogy  between  the  stability  of  rotating 
flows  and  the  stability  of  a  stratincd  flukl  at  rest  in  a  gravitational  field,  the  analogue  I'f 
<?  in  fact  being  the  square  of  the  Brunt-Vaisaila  frequency. 

In  subsequent  sections  a  dynamical  system  is  established  to  discuss  the  generation  of 
vertical  cells  in.sidc  vortices.  The  sea  water  is  assumed  t<'  be  Boussinesq  anil  inviscid. 
Therefore,  the  system  excludes  the  damping  ciTcci  of  frictional  forces  on  the  vortex. 
Furthermore,  for  the  sake  of  simplicity,  we  also  assume  no  vertical  vcUk'Iiv  at  the  top  and 
the  bottom  of  the  vortex  and  no  interaction  between  the  vortex  and  the  ambient  flow. 


2.  dynamical  system 

2.1  Total  Row 

A  cylindrical  coordinaic  system  (r, z)  is  used  in  this  study  with  the  rotating  axis  of 
the  vortex  as  the  center  of  the  sj-stem  (Fig. I).  The  vertical  cinirdinaic  is  positive  upward 
and  zero  at  ihc  bottom  of  the  vortex. 
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Fig.  I.  The  cylindrical  c<M>rdinaics  (r.  A, a)  with  the  vertical  axis  at  the  center  of  the 
vortex. 

The  Bou.ssincsq  approximation  is  a.ssumcd  for  our  system,  is  the  characteristic 
density  of  the  sea  water.  The  total  density  and  pressure  can  hr  decomposed  into 

/'  =  f>(i  +  (>-  (2/7) 

P‘ 

where  p,  is  the  air  pressure  at  the  ivccan  surface,  p.  is  the  dynamieal  pressure,  anil  p.  is  the 
density  deviation  from  the  characteristic  value.  Under  the  assumpliim  that  ihc  moiion 
is  axisymmctric  and  inviscid.  the  governing  cqu.ations  arc  written  in  the  eyiindrical  coor¬ 
dinates  as 


r 


)r  = 


I  PP’ 
f'o  Or 


(:>«) 


(3/0 
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dw  _  _  _ I 

dt  ~  Po  ^  Pn  ,?r 


(3/-) 


^  c{rw)  _ 
Pr  '  (iz  ' 


(3//) 


dp. 

-^  =  G  (3^) 

where  w.  i'.  w  arc  radial,  tangencial.  an<l  vertical  componenttortlic  flow  ncld.  rcwpcciivcly. 
Q  is  the  thermally  induced  rate  of  change  of  density.  However,  the  dynamically  driven 
vertical  cells  and  their  contribution  to  the  ilccp  convection  are  the  main  items;  ihcreft're. 
Q  Is  set  to  be  zero  in  this  research.  The  operator  djdt  is  defined  by 

d  ^  ^  ,  '3  « 


2.2  Basic  Flow 

A  steady  axLsymmetric  and  gr.idicntly  balanced  circular  vortex  with  tieprh  H  and  ra¬ 
dius  R.  is  taken  to  be  the  basic  flow,  i.e., 

i7  =  0.  v  =  ffr.  ^).  v?t=0,  for0<r</?,  0  <  r  < /y  fS) 

The  momentum  equations  are 


I  ^P 
Pn  or 


(fw) 


(6/d 


where  V  is  the  mean  tangential  current,  /  is  the  Coriolis  paramcicr.  The  differentiation 
of  (6a)  with  respect  to  z  minus  the  slifferentiaiion  of  (6b)  with  respect  to  r.  leads  to  the 
thermal  wind  relation  for  the  basic  flow 


S  dp 
Pn  dr 


(7<r) 


Therefore,  the  basic  flow  is  the  gradient  balanced  flow.  Including  the  earth  rotation,  the 
angular  velocity  of  water  parcels  at  a  radius  of  r  Ls 
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The  square  of  the  angular  momentum  is 


Af  =  (nrV  =  (X.  -L,v 

The  thermal  wind  rciaiion  (7a)  becomes 

1  <11/  _  y  ^ 

re  Po  Pr 


Oh) 


2.3  Disturbances 

Suppose  that  axisymmetrk  disturbances  arc  generated  from  ihe  basic  fiovt.Thc  total 
flow  fickl  is  then  a  combination  of  the  iwo  parts;  basic  fiow  an<l  tlisttirbnnccs, 

=  w',  r  =  !•  +  r'.  w  =  w'.  =  'p  •*■(>'.  f>.  —  /i  ^  p'  (8) 

Substituting  (S)  into  (3a-c)  and  taking  into  account  ((ia.b).  and  assuming  small  ampli¬ 
tudes  of  disturbances,  a  set  of  equations  for  the  disturbances  is  then  obtained 


t  ^  Po  Pr 

Oc) 

'  ct  cz 

Oh) 

fiw' _ ^  1  hp' 

Pi  P(\  ^  Pa  rr 

(9r) 

P{ru‘)  <)(w') 

-  a  ~  ^ 

tir  cz 

Od) 

cp'  rifi  .(fp, 

~ - r  1/  -  W  =  0 

rl  rr  cz 

Oc) 

The  venical  circulation  generated  inside  the  vortex  can  be  depicted  by  the 
streamfunction 
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t/  “ 


I  0<J' 


w  ~ 


Ih 


(10) 


Eliminstiiig  v'.p'.p'  from  (9.i-c).  we  uhiain  an  equation  for  the  s!tcnmruncii«>n  >]i. 


— i  ■  I  1  ,  S  ,  A  B  , 

Sr-  ^  ^  ^  ,1,2  J  r  pr  "  '  cz  ^ 


C 


f?«// 

»)? 


B  <V 

-r^)  =  « 

cr 


(ID 


which  is  ihc  same  equation  used  by  Ooyama  (1966)  and  Charncy  and  where 


S  >2p 
^'’n  cz 


(I2fl) 


I  _2v_  .  _ 1?  t'/' 

’  ^  (?r  ~  ^11  (?r 


(12/.) 


Cs 


(I2<') 


represent  three  clifTcrcnt  frequencies:  ^  A  is  the  buoyancy  frequency,  ^  B  is  the  baroclinic 
Irequency.  and  Ls  the  gcneralizctl  inertial  frequency.  Here,  the  mean  density  p  is  a 
(unction  of  r  and  :  From  (12a)  and  ( 12b)  we  can  get  the  relationship  between  -1  and  B: 


dA  __  dB  _  ]  B'  u 
dr  dz  p,' 


(12./) 


If  the  oce.an  has  constant  stratification  out.sklc  the  sortex.  (he  integration  ol'(l2d)  with 
respect  to  r  leads  to 


.1  =  A"  +  D 


( 1  ?« ( 


D  = 


2 

r 


(fn* 


(l3/>) 


where  A'  is  the  characteristic  value  of  the  Brunt-Vai.sala  frequency,  which  is  a  constant  in 
this  research.  The  time  scale  for  deep  convection  is  generally  much  longer  than  A  '  and 


/ i-c.. 
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I 


di 


The  lime  rate  or  change  shoiili!  be  neglected  against  both  the  buoyancy  Trcqucncy  and  the 
inertial  frequency.  Therefore,  the  equation  for  sircamfunction  (  I  I)  is  simplified  into 


A^!L)  cl  c  ,  B 

,7r  '  ''  nr  ~  nz  ^  rz  "  r7r 


(14) 


which  shows  the  vertical  ceils  in  (c.  z)  section  generated  by  ihe  combined  effect  of 
buoyancy,  baroclinicity.  and  inertia  of  a  vortex,  which  is  called  the  symmetric  instability. 
Here,  we  .arc  interested  in  the  intrinsic  physical  mechanisms  generating  \erncal  cells  inside 
a  vortex.  Interaction  between  the  vortex  and  the  ambient  fltnv  is  beyond  the  scope  of  this 
research-  Therefore,  it  is  reasonable  to  assume  that  n<i  flow  can  across  the  ctlgc  of  the 
vortex.  The  radial  velocity  at  the  center  of  the  voncx  should  be  ?cro.  Furthermore,  we 
also  assume  no  vertical  velocity  at  the  top  and  the  btutom  of  ihr  voncx  for  the  sake  of 
simplicity.  Thus,  the  boundary  conditions  for  (14)  arc 

1^  =  0.  at  r  =  0.  r=R.  ?  =  0.  (15) 


3.  NONDIMENSIONALIZATION 

If  the  independent  variables  (r.  zl.  mean  tangential  vchKity.  .squares  of  the  three  fre¬ 
quencies  (P.  C,  D),  and  the  streamfunction  arc  nondimcnslonalizcd  by 


(r.  r.  v)  =  IRr.  Hi,  Kv), 


( B,  C .  D)  =  f(Ro  B.  C.  Ro  h).  <l> 
H  /y- 


I  H-iU  (Kv) 


where  I'  is  the  characteristic  tangential  vchxrity  of  the  vortex,  and 


«  =  -  (1  -f 

r  CZ 

(\7c) 

C  =  (l  +  /(o  ^  4- -2^  ) 

r  r  cr 

(nh) 

D  =  -  f'[(l  4 

J-  r  ffz'^  r  CZ 

(17c) 
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where  Rn  =  i^l/R  is  the  Rnsshy  number.  The  nondimcnsional  form  ut"  the  equation  (14) 
is 


(5.-/  -  RoD] 


(7^ 

^7- 


-  d^ij/ 
I  2RoB- 


^ri^z 


—  I  iiiii 
-{Bu  +  RoP)-l:r  — 
r  Br 


cC 


SB 


fi^i-Ro 
Sz  (ir 


a  Si!/ 

-  R>'-^)-^  =  o 

(•  f?’ 


(IS) 


where  Bit  =  (NHIfRY.  is  the  Burger  number,  which  is  the  square  of  the  ratio  of  the 
Rossby  radius  »>f  rleformaiion  to  the  radius  of  the  vortex.  Ft*r  .1  given  vortex,  the  m.ag- 
nitude  of  Bit  indie.ates  the  s! ratification.  The  smaller  the  Bu.  the  weaker  the  stratification. 
The  nonilimensional  form  of  the  bounil.ary  conditions  (l.S’)  is 

4/^0,  ai7-n,  ?=l.  ?=0.  ?=l  (]<)} 


4.  VERTICAL  CELLS  DRIVEN  BY  A  GAUSSIAN  VORTE.X 

A  cyclonic  vortex  with  Gaussian  distribution  in  both  radial  an<l  \ crlic.il 

v(r,c)  =cxp[-/7i,(r  -  -ip)’]  (20) 

is  chosen  as  a  basic  flow  in  this  study,  where  m.,  tn.  arc  parameters  intlicating  the  strength 
of  the  shear  in  radiai  and  vertical  directions.  Equation  (IS)  with  the  bmindary  condition 
(l‘>)  and  mean  Bow  (20)  is  integrated  numerically.  The  successive  over-relaxation  method 
is  used  for  the  integration.  The  solutions  show  the  generation  t»f  vertical  celis  inside  the 
vortex  even  in  the  relatively  strong  stratification  case  (Bu  -  I).  The  other  three  movlcl 
parameters  Rn.  m,.  m,  arc  fotind  to  be  important  in  affecting  the  structure  of  the  vertical 
cells  insirlc  the  vortex.  For  a  given  vortex  size,  each  parameter,  taking  only  ivvo  different 
values,  represents  the  two  different  physical  conditions,  such  as.  R/t  —  0.1.  for  a  relatively 
weak  vortex,  and  Ro  =  0.2.  for  a  relatively  strong  vtvrtcx:  m.  =  0.  for  the  absence  of  ratlial 
shear,  and  m,  =  4.  for  a  rciaiivcly  strong  r.idial  shear;  an<l  m,  —  0.  for  the  absence  of  ver¬ 
tical  shear,  anil  m,  =  4.  for  a  relatively  strong  vertical  shear. 


Fic.2  show  ihc  sircamrunction  in  a  radial-vertical  section  To:  Bu  1  and  Ro  -0.1: 

(a)  m,  =  0.  m,  =  0,  (/?)  m,  =  4.  m,  =  0.  (r)  m.  =0.  m,=  4.  (r/j  m,  =  4.  rri,  =  4.  These 
eases  refer  to  a  rclaiively  strong  siraiification  and  a  relatively  weak  vortex. 

Several  rcsulus  can  be  easily  drawn  from  Fig.2: 

fa)  The  vertical  cells  <lo  exist  inside  the  vortex  for  all  eases,  even  m  the  ease  without 
radial  and  vertical  shears  of  Ihc  mean  tangential  velocity  as  shown  in  Fig.23.  This 
stronglv  suggests  that  the  square  of  the  generalized  inertial  ricqucncy.  r(r).  is  a  funda¬ 
mental  factor  in  generating  vertical  circulation  inside  the  v<irtcx  (since  B  =  0  and  D  -  0). 

(b)  Comparison  of  Fig.2a  with  Fip.2b  shows  that  without  the  vertical  shear  on  the 
mean  t.angcntial  velocity,  the  cfrcct  of  the  radial  shcai  (m,  ^  0)  on  the  vcrtic.'il  cells  is  very 
small. 

(c)  Comparison  of  Fip.2b  with  Fig.2d  indicates  thm  with  vcntcal  shear  on  the  mean 
tangential  velocity,  the  effect  of  Ihc  radial  shear  on  the  vetiical  cells  becomes  evident. 
The  upper  cell  penetrates  it'  the  base  of  the  vortex,  and  the  tvvi>-ccll  structure  becomes  a 
ihrcc-ccll  structure. 

fd)  Comparison  of  Fig.2c  with  Fig.2a  leads  to  the  result  that  the  baroclinity  always 
inicnsirics  the  strength  of  the  vertical  cells. 

Fig.3  shows  the  strcamfunciion  in  a  radial-vertical  .section  for  Bn  =  1  and  Rn  =  0.2; 
(a)  m,  =  0.  m,  =  0.  (f>)  m.  =  4.  =  0.  (e)  w,  =  0,  m,  =  4,  (rf)  m.  =  4.  m.  =  4.  These 

ca.scs  refer  to  a  relatively  strong  siraliftcation  and  a  relatively  str<*ng  vortex.  Comparison 
of  Fig..'  with  Fig.2  shows  that  the  increase  of  the  Rossbv  nuinlvcr  will  generally 
strengthen  the  vertical  cells.  Therefore,  for  a  given  size  of  the  vortex,  the  faster  the  rota¬ 
tion  of  the  vortex,  the  stronger  the  vertical  cells  inside  the  vortex. 


Fig.  2.  The  sircam fund  ion  in  ratlial-vcriical  section  for  Sm=I  anti  /?/)  =  ().! 
(a)  m,  =  0.  m.  =  0.  (A)  m.  =  4,  /n,  =  0.  (c)  m.  =  0.  w.—  4.  (</)  m,  =  4,  w.  =  4. 


5.  TIME  RATE  OF  DENSITY  INSIDE  THE  VORTEX 


Define  a  nonditncnsional  density  pcnurbaiion  averaged  »'\cr  the  horizontal  plane  in¬ 
side  the  vortex  as  follows: 


n 


rdr^ 


•■(1 


(21) 


Time  t  is  nondimcnsionalizcd  by 


f=Cr')— (22) 

2Fr-/> 

where  Fr  =  is  the  Froude  number,  and  =  R  //  ,  is  the  aspect  ratio.  Inlcgr-’iiing 

the  density  equation  (9c)  over  the  horizontal  plane  within  the  vortex,  anti  using  the  non- 
dimensional  forms  of  parameters  (l7a.c),  a  nondimcnsional  equation  tif  time  rate  of 
change  of  htirizontaliy  averaged  dcn.sity  is  obtained 


a<T 

n7 


(2.^1 


The  derivative,  dalat.  is  computed  for  the  ca.se:  Bii=\.  /?o  —  0.1.  //;,  = '4.  w.  =  I 
(Fig. 2d).  The  effect  of  the  vortical  motion  on  the  density  redistribution  inside  the  vortex 
Ls  shown  in  Fig.4,  where  the  density  decreases  in  the  lower  part  {}  <  ().2~)  of  the  vortex, 
and  increases  in  the  upper  part  (e>0.27)  of  the  vortex.  Tltis  pioccss  rctiuces  the 
stratification  and  tends  to  make  the  ilcnsity  vertically  uniform  inside  the  vortex.  We  view 
this  priKcss  as  one  contribution  to  the  precomlitioning  of  the  open  ocean  deep  eonvcctiim. 
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6.  CONCLUSIONS 

(a)  This  research  is  intended  to  bring  to  the  attention  of  the  (Kcanographic  community 
the  importance  of  symmetric  instability  in  vortex  dynamics,  and  in  generating  convection. 
Solutions  show  generation  of  a  variety  of  convective  cells  inside  the  \ortcx  by  this  mech¬ 
anism. 

(b)  The  vertical  cells  can  be  generated  in  the  vortex  with  relatively  strong  S[ratificatK)n 
and  without  any  radial  or  vertical  shear.  The  rotation  itself  svill  intiucc  the  vertical  cir¬ 
culation. 

(cl  The  baroclinicity  always  intcnsirics  the  vertical  circulation.  However,  the  clTcct 
of  the  radial  shear  on  the  \cnical  cells  Isccomcs  import.nnt  only  when  it  is  associated  with 
vertical  shear. 

(vl)  The  stronger  the  rotation  of  the  vortex,  the  stronger  the  vertical  cells  inside  the 
vortex. 

(c)  The  lime  scale  for  the  density  redistribution  by  this  dynamically  driven  overturn¬ 
ing  IS  inversely  proportional  to  the  s<|uarc  of  the  Froude  number.  For  numerous  olsscr- 
vaiions  of  vortices  in  the  ocean,  this  time  scale  is  several  weeks.  This  scale  is  coincident 
with  the  time  period  needed  for  the  preconditioning. 

(0  The  .solution  here  refers  to  the  Gaussian  type  cyclonic  vortex.  For  other  kinds  of 
vortices,  the  solution  may  vary. 

(g)  The  solution  is  obtained  from  assuming  rigid  boundaries  Poth  .it  the  top  and  the 
Nstiom  of  the  vortex.  If  free  boundary  conditions  or  other  kintls  of  boundary  conditions 
arc  usctl.  the  shape  of  the  v  ertical  cells  will  change,  especially  near  the  boundaries,  but  the 
basic  physics  will  remain. 
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Galeti 


Fig.  4.  Vertical  dependency  of  for  the  vortex  with  Bu=\.  Ro  =  0.\.  w,  =4.  and 

m,  —  A. 

The  lime  scale  for  the  density  redistribution  process  inside  the  vortex  can  he  roughly 

estimated  as  follows.  If  the  tangential  velocity  for  the  vortex  f  is  5fi  cm-s.  and  if  the 

thickness  of  the  vortex  H  is  on  the  order  of  I  km.  the  square  of  the  Froude  number 

will  he  0.2.‘<  x  lO  ^  The  aspect  ratio  fi  is  on  the  order  of  0.1,  and  the  change  of  the  non- 

dimcnsi<*nal  density  pcrtiirhation  A"’  is  on  the  order  of  10  and  the  integration. 

/=  \'(Bfiti/ldz  ~  b^^!vr)dr  ,  is  on  the  order  of  .*>  (Fig.  41.  Therefore,  the  lime  scale  for 
•■o 

this  density  redistribution  process  Ls 

7'-'/"’ - ¥- - -40r'.  (24) 

2Fr-!>0{l) 

which  is  on  the  order  of  few  weeks.  This  scale  is  coincident  with  the  time  period  (few 
weeks)  needed  for  the  preconditioning  which  creates  a  region  of  very  weak  .static  .stability 
within  the  voncx  (Kill worth.  1985). 
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